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We investigate transition of a one-dimensional interacting Bose gas from a strongly repulsive 
regime to a strongly attractive regime, where a stable highly excited state known as the super 
Tonks-Girardeau gas was experimentally realized very recently. By solving exact dynamics of the 
integrable Lieb-Liniger Bose gas, we demonstrate that such an excited gas state can be a very 
stable dynamic state. Furthermore we calculate the breathing mode of the super Tonks-Girardeau 
gas which is found to be in good agreement with experimental observation. Our results show that 
the highly excited super Tonks-Girardeau gas phase can be well understood from the fundamental 
theory of the solvable Bose gas. 
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Introduction. — Recent experimental advances with ul- 
tracold atoms have opened up a new avenue for the study 
of one-dimensional (ID) strongly correlated many-body 
systems which continue to inspire significant develop- 
ments in physics [lrQ. Exquisite tunability of ID effec- 
tive interacting strength via Feshbach resonance or the 
confinement-induced resonance 043 allows the simula- 
tion of low-dimensional quantum gases with arbitrary 
interaction strength. The experimental realization of 
Tonks-Girardeau (TG) gases p§, [Ij has provided a direct 
test on the fermionization of Bose gas in the strongly in- 
teracting limit [1, ij . Very recently, Haller et al [l(| have 
made a new experimental breakthrough on the realiza- 
tion of the stable highly excited gas-like phase called the 
super-Tonks-Girardeau (STG) gas, which was predicted 
theoretically and verified from the integrable inter- 
acting Bose gas with attractive interactions [l2[ • Their 
experimental results have stimulated further theoretical 
study of such novel states with strong correlations in ID 
many-body systems which was also found theoret- 
ical correspondence in a strongly attractive Fermi gas 
without polarization [l4j . 

In general, experimental realization of a stable excited 
state is difficult and a little counter-intuition since it has 
no analogue in traditional solid state systems, where a 
pure excited state is not stable due to the energy dis- 
sipation between the system and the environment. The 
experimental progress [lfj| opens the door for searching 
novel quantum states in quantum gases and studying 
physical properties of stable excited phases, which pro- 
vides a promising new area of activity in cold atoms. 
However, the metastability, dynamics and statistical sig- 
nature of the excited states are still far from under- 
stood. Current theoretical understanding of the STG 
gas is based on intuitive explanation that the STG gas 
inherits hard core behavior from the repulsive TG gas 
which prevents the collapse to its true cluster ground 
state (GS) [TTI - [T7l |- As a matter of fact, the STG gas is 



obtained from the TG gas by a sudden switch of interac- 
tions, this process is in principle a dynamics problem for 
the interacting boson model. Therefore a comprehensive 
understanding of the STG gas is highly desirable and it 
is also important to show that such a stable state really 
follows from the fundamental quantum dynamics of an 
exactly solvable model. 

In this work, we demonstrate how the GS of the repul- 
sive TG gas translates to the highly excited state of the 
strongly attractive TG gas through a sudden switch of 
interactions. We also show that the highly excited STG 
phase is stable and the transition rate to the GS is almost 
completely suppressed in the strongly interacting regime. 
Furthermore, we calculate the energy of the STG gas by 
solving the exact Bethe ansatz equations (BAEs) and de- 
termine the collective excitation modes of the STG gas 
with the help of the local density approximation (LDA), 
which agrees quantitatively with the experimental data. 

Interacting boson model. — We consider a quasi-lD sys- 
tem with TV bosons tightly confined in an elongated trap 
which is described by an effective ID Hamiltonian 

N h 2 d 2 

i—l L i<j 

with g = —2h 2 /(maid) being the effective ID interac- 
tion strength and aid the effective ID scattering length 
0, H[ . As a simple but fundamental integrable interact- 
ing boson model, the exact results for the Lieb-Liniger 
(LL) model |T]) play an important role in understanding 
quantum statistical and many-body correlation effects of 
a ID bosonic quantum gas |9l. [l8l - l21f . The model (H) is 
exactly solved by the Bethe ansatz method with eigen- 
functions taken as superpositions of plane waves over all 
permutations of the momenta 
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can be determined by solving the BAEs 
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For the solutions with Ij = (N + l)/2—i, the correspond- 
ing state is just the STG gas state [T^| . In the strongly 
attractive interaction limit (c — > -co), the quasi- 
momentum distribution of real roots determined from ([5]) 
reduces to the free Fermi distribution with a minimum 
separation between two quasimomenta 2ir/L. As the at- 
traction strength decreases, the gas-like highly excited 
state has a stronger anyonic-like pressure than the Fermi 
pressure. For strongly attractive interaction (c — > — oo), 
the energy of the STG gas state [l2| is given by E$tg = 
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FIG. 1: (color online) Quasi-momentum distributions for the 
ground state of the repulsive and the STG gas phase of the 
attractive Bose gas with different values of 7. 



continuously connected to the GS energy of the TG 
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as \c\ -> 00 sj. 

Quench dynamics — Suppose that the initial state 
\^ in (X,t — 0)) = \ip (X,c)) is prepared in the strongly 
repulsive limit c> 1, after a quick switch from the repul- 
sive regime to the attractive regime with d < 0, the wave- 
are functions of the two- function \V(X, t)) = e im \^ t (X, c)) can be calculated via 



where the wave functions are defined in the domain 
x\ < X2 < ■ ■ ■ < xn and the coefficients Ap = 



particle scattering phase shifts with (— I ) p = ±1 for odd 
or even P. The quasi-momenta kj are determined by the 
BAEs [|] 
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Here the coupling constant g = h 2 c/m with interaction 
strength c — determined by the effective ID scat- 

tering length aid- The total momenta of the system is 



given by K — 53j=i % an< ^ the eigenenergies are given 

byjB = |lvf , kl 

For repulsive interaction c > 0, all solutions kj of © 
with j — 1, • • • , A are real numbers. Alternatively the 
solutions follow by solving the log form 
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of the BAEs. The GS solution corresponds to Ij — (N + 
l)/2 — j. On the other hand, for attractive interaction 
c < 0, the GS solution to the BAEs follows from the N- 
string solution kj = Xj + i(N + 1 — 2j)(c/2 + 5j), where 
j = 1 , ■ • • , N and 5j are small deviations which fall off 
exponentially to zero as c — > — 00. Correspondingly, the 
GS is described by a cluster state [l5[ with the energy 
E = -^c 2 A(A 2 - 1). 

STG phase — Although the GS solution corresponds 
to the A-string solution, the BAEs ((3]) still have real 
solutions even for c < 0, which however correspond to 
some highly excited state of the attractive Bose gas and 



where c„ = (ip n (X,d)\i/)o(X,c)) with ip n (X,d) repre- 
senting the n-th eigenstate of the LL model with at- 
tractive interaction strength d and X the abbrevia- 
tion of xi,- ■ ■ ,xjy. From ([6]), it is straightforward that 
the probability for the state after quench staying in a 
STG phase is given by \{4>stg(X, c')\M x , c))| 2 , which 
also represents the transition probability from an ini- 
tial TG phase to the final STG phase. The wavefunc- 
tion tpsTG is determined by the Bethe ansatz wave func- 
tion @ with kj determined by the solutions of (J5J) with 
Ij = (N + l)/2 — j. Similarly, the transition probability 
from the initial TG state to the cluster state is given by 
\ (ip duster (X, d)\ipo(X, c))| 2 . Since the Bethe ansatz wave 
functions tp n (X, c) can be exactly determined by solving 
BAEs, in principle we can calculate the transition prob- 
abilities exactly. 

For simplicity, we first consider the case with d = — c, 
i.e., the system is initially in the GS of the LL gas with 
c > 0, and then suddenly switch to the attractive regime 
with strength — c. We note that the BAE solutions for 
the repulsive LL gas and that for the STG gas corre- 
spond to the same set of Ij according to Eqs. (d]) and 
(O. In the limit c = —d — > 00, the solutions are given 
by kj — Ij2ir/L, i.e. they are exactly the same. In 
Fig. I, for particle number N — 9, we show the BAE 
solutions for the repulsive LL gas and the STG gas with 
different values of 7 = c/p. In the strongly interact- 
ing regime, i.e. |— y | 3> 1, the quasi-momentum distri- 
butions for the repulsive LL gas and the STG gas ap- 
proach the free fermion "orbitals" from different sides. 
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FIG. 2: (color online) Transition probabilities from the TG 

gas to STG phase. Inset: the transition probabilities from FIG. 3: e(-y) and 52(7) versus | — / 1 . 

TG gas to the cluster state. 



As I7I decreases, the deviation from the free-fermion dis- 
tribution become more obvious. From the change of 
momentum distribution, we can expect that the overlap 
between the GS wavefunction of the strongly repulsive 
LL gas and the eigen- wavefunction for a STG gas, i.e., 
\{i/}STG(X,-c)\ipo(X,c))\, is close to 1 for very large | 1 
and exactly 1 when I7I — > 00. 

To see how the transition rate changes with respect 
to 7, we calculate the overlap of the wavefunctions be- 
tween the repulsive and attractive regimes. In order to 
conceive the signature of the overlap between the two 
regimes, we consider a system with N = 4 and calculate 
the transition rate from the initial repulsive LL gas to the 
STG gas with different In Fig. 2, we show the tran- 
sition rate from an initial repulsive gas with 7 = 200, 
to a final STG phase with different values 7 < 0. As 
expected, the transition probability to the STG phase is 
close to 1 after switching to the strongly attractive regime 
whereas the probability for dynamically falling into the 
cluster state is almost completely suppressed. However, 
when the system is switched into the weakly attractive 
regime, the STG phase is no longer stable and the prob- 
ability for falling into the cluster state increases quickly. 
The transition rate for the larger system does not change 
qualitatively, however the calculation for a large system 
is a very time-consuming task due to the calculation of 
multidimensional integrals. Our calculations give a clear 
signature of the metastable STG gas against collapsing 
into the cluster state. 

Properties for the STG gas — As shown in the above 
calculation, the transition rate to the STG phase is very 
high in the strongly interacting regime. Therefore such 
a highly excited state with a stronger pressure than the 
Fermi pressure could be reached through switching the 
interaction from strong repulsion into strong attraction. 
We note that such pressure is inherited from the repulsive 
gas, i.e. transmutes the statistical kinetic energy into the 
attractive STG phase. In the thermodynamic limit with 
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FIG. 4: (color online) Breathing mode of the STG gas. The 
experiment data with error bars and Monte Carlo result are 
reproduced from Fig. 3a of Ref. [T(j. 



N, L — > 00, the eigen energy of the STG phase can be rep- 
resented in the integral form: E S tg/L = J k 2 p(k)^ 



and p(k) 
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gration limit Q is determined by p 
making replacement of variables k 

g(x) — p{Qx) according to 0, the energy per atom can 
be expressed as e (p) = j^p 2 e (7), which can be obtained 
by solving the system of equations 
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and A = I7I f -. g(x)dx. The above equations are 
very similar to Lieb-Liniger's solutions for the repul- 
sive boson gas except the minus sign in Eq. ([5]). By 
numerically solving the integral equation, we can ob- 
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tain e(7) and thus the energy density e (p) and the 
chemical potential p,(p) = d p [pe(p)}. For the STG 
phase in the strongly attractive regime, we find e(7) ~ 

In Fig. 3, we show e("f) for different values of 7, which 
is found to be well fitted in the whole attractive regime 
by a rational function e (7) = ^§-(1 + Pi I7I + P27 2 + 

P | 7 | 3 /4)/(l + gi | 7 | + 927 2 + P |7| 3 ) with pi = 0.075, 
p 2 = 0.013, gi = 0.227, q 2 = 0.034, p = 0.004. We also 
calculate the local two-particle correlation function 172 = 
(\]/t\j/t\j/\]/) which can be obtained by 32(7) = p 2 de^) / dj 
[20| . The STG phase exhibits stronger local correlation 
than the repulsive Bose gas. For the STG gas in a har- 
monic trap with V ex t = mui 2 x 2 /2, we can determine the 
density distribution of the STG gas using LDA pll. I22T] . 

According to the LDA, one can assume that the system 
is in local equilibrium at each point x in the external trap. 
The density distribution of the STG gas is then obtained 
from the local equation of state p[p(x)] = pa — V ex t{x) un- 
der Thomas- Fermi approximation with the normalization 
condition J dxp(x) — N. Using the density distribution, 
one can then calculate the mean square radius of the 
trapped STG gas via (x 2 ) — J p(x)x 2 dx/N. Following 
Ref . [ll|, Hll , we can calculate the frequency of the low- 
est breathing mode given by w 2 = -2(x 2 )/(d(x 2 )/du 2 x ). 
Comparing with the Monte Carlo simulations [111 ], our re- 
sult (solid line) based on the exact solutions of BAEs has 
a pronounced peak with a maximum of breathing mode 
frequency w 2 /u> 2 reaching the value about 4.5 which co- 
incides with the experimental data, see Fig |4j Here the 
discrepancy between the theoretical prediction from the 
LL model and the experimental data inspires us to search 
for a new theoretical framework capable of understand- 
ing the quantum dynamics of the excited quantum gas 
phase. 

To this goal, we notice that the resonance regime of 
switching the interaction from strong repulsion to strong 
attraction can be identified as anyonic type of interaction 
c/ cos(k/2) [23|,[24j]. We also understand that strongly in- 



teracting bosons have lost their bosonic signature due 
to the transmutation between statistical and dynami- 
cal interactions [3]. This naturally suggests us to con- 
sider the integrable model of anyons as a fitting model 
to the experiment on the STG gas. The dynamic pro- 
cess of switching the interaction from strong repulsion to 
strong attraction may be viewed as the change of quan- 
tum statistics from less to more exclusive than Fermi 
statistics. In order to conceive this nature of the STG 
gas, we also calculate the breathing mode of the inte- 
grable model of anyons with n — 3.93 greater than the 
free Fermi case k — it, see dash-dotted line. The good 
agreement with experiment data implies that the any- 
onic signature may relate to the statistical signature of 
the highly excited state of interacting bosons. However, 
we notice that in the weak coupling regime the involve- 
ment of cluster states of few particles [l3| is the main 
reason for the discrepancy between theoretical prediction 
and experimental data. When the system is not in the 
strongly attractive regime, the STG phase is no longer 
stable, see Fig|4] 

Summary — In summary, we have investigated the 
transition from the GS of strongly repulsive bosons to the 
gas-like highly excited state of attractive bosons through 
a switch of interaction. By solving the quench dynamics 
problem of the interacting bosons, we have shown that 
the gas-like excited state is stable in the strongly attrac- 
tive regime. Our exact result for the frequency of the low- 
est breathing mode obtained from the BAEs for interact- 
ing bosons and anyons are in reasonably good agreement 
with experimental observation. This suggests a statisti- 
cal signature of the STG gas which may be viewed as an 
anyon-like gas. 
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